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System and Method for Analyzing and Imaging 
Three-dimensional Volume Data Sets 



Background of the Invention 



Field of the Invention 

Tlie present invention relates generally to imaging of three-dimensional ("3D") 
volume data sets. More particulaiiy , the present invention relates to fast visualization and 
analysis of structures within 3D volume data sets. 

Related Art 

Many fields of endeavor require the analysis and imaging of three-dimensionial 
("3D") volume data sets. For example, in the medical field, a CAT (computerized axial 
tomography) scanner or a magnetic resonance imaging (MRI) device is used to produce 
a "picture" or diagnostic image of some part of a patient* s body. The scaimer or MRI 
device generates a 3D volume data set that needs to be imaged or displayed so that 
medical personnel can analyze the image and form a diagnosis. 

Three-dimensional volume data sets are also iised in various fields of endeavor 
relating to the earth sciences. Seismic soimding is one method for exploring the 
subsurfiu:e geology of flie earth. An imderground explosion or earthquake excites seismic 
waves, similar to low fi^quency sound waves, that travel below the suifece of earth and 
are detected by seismogr^hs. The seismographs record the time of arrival of the seismic 
waves, both du-ect and reflected waves. Knowing the time and place of the explosion or 
earthquake, the time of travel of the waves through the interior can be calculated and used 
to measure the velocity of the waves in the interior. A shnilar technique can be used for 
offshore oil and gas exploration. In ofiEshore exploration, a ship tows a sound source and 
underwater hydrophones. Low fi-equency (e.g., 50 Hz) soimd waves are generated by, 
for example, a pneumatic device that works like a balloon burst. The soimds bounce off 
rock layers below the sea floor and are picked up by the hydrophones. In Ihis manner, 
subsurface sedimentary structures that trap oil, such as faults, folds, and domes, are 
"mapped*' by the reflected waves. The data is processed to produce 3D volume data sets 
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that include a reflection or seismic amplitude datavalue at specified (x, y, z) locations 
within a geographic space. 

A 3D volume data set is made up of "voxels" or volume elements. Each voxel 
has a numeric value for some measured or calculated property, e.g., seismic amplitude 
of the volume at that location. One conventional approach to generating an image of a 
3D volume data set is to cross-section the 3D volume data set into a plurality of two- 
dimensional ("2D") cross-sections or slices. The image of the 3D volume data set is then 
buih as a composite of the 2D slices. For example, the image of the 3D volume data set 
is generated by staclcing the 2D slices in order, back-to-front, and then composited into 
a complete image. The user sees the image being built layer by layer as the composite 
grows. Although the user can see the internal organization or structure of the volume as 
the composite unage grows, the traditional slice and composite technique is typically 
slow, particularly when very large 3D volume data sets are being used. Additionally, the 
sUce and composite technique clutters the user's field of view with ex-traneous 
information, and interferes with the user's ability to accurately visualize and interpret 
features inherent in tlie 3D volume data set 

Computer software has been developed specifically for imaging 3D seismic data 
sets for the oil and gas industry. Examples of such conventional computer programs 
include VoxelGeo, available from Paradigm Geophysical, Houston, TX, Seis Works and 
EarthCube, available fi-om Landmark Graphics Corporation, and lESX, available fi'om 
GeoQuest Such conventional computer programs have numerous deficiencies that 
precludeauserfiromquickly and accurately visualizing and inteq^ 
in a 3D seismic data set Conventional computer programs for visualizing and 
interpretmg 3D seismic data operate on the fiiU 3D volume of seismic data. 
Consequently, every time a change is made, such as a change to the transparency or 
opacity settings, the fiill 3D volume of seismic data must be processed, and the image 
re-drawn. Even when such programs are run on highly efficient graphics 
supercomputers, the delay or lag in re-drawing the image is perceptible to the user. For 
a 3D volume containing 500 megabytes of seismic data, it can take on the order of 30-45 
seconds for conventional programs to re-draw the complete image (fi:ame rate of 0.03 to 
0.02 fi:ames per second, respectively). During tiie 30-45 second delay time, the mind of 
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the user loses focus on the feature of interest, making it difficult to completely and 
properly analyze tlie seismic data. 

Some conventional 3D seismic interpretation programs provide the capability to 
visualize and interpret a piece of the full 3D volume of seismic data. The user identifies 
the coordinates of the selected piece via a menu command. An image of the selected 
piece is drawn. The selected piece can then be rotated, if desired, at that location. 
However, to look at a different piece of the full 3D volume of seismic data, such as to 
follow a geologic feature that has been tentatively identified, the image must be 
interrupted, a new location or coordinates for the different piece is entered, and a new 
image is drawn containing the different piece. The interruption in the displayed image 
makes it difficult for the user to visualize any continuity between the two pieces of the 
fiill 3D volume of seismic data that have been imaged. This impedes the user's ability 
to interpret and identify the geologic features that are inherent in the full 3D volume of 
seismic data. Additionally, even though only a piece of the full 3D volume of seismic 
data is being visibly displayed, conventional 3D seismic interpretation programs continue 
processing the full 3D volxmie of seismic data to draw the image, thereby slowing the 
display of the image to the user. 

Conventional 3D seismic interpretation programs provide the capability to "auto 
pick" and identify points that satisfy a voxel selection algorithm. However, these 
programs typically iterate through the full 3D volume of seismic data to identify the 
points that satisfy tiie voxel selection algorithm. This is time consuming even on a high 
speed graphics si5)ercomputer. Additionally, conventional 3D seismic interpretation 
programs do not provide the capability to directly delete from Ae collection of picked 
voxels. Hie only way to "eliminate" points from the collection of picked voxels using 
conventional 3D seismic interpretation programs is to repeatedly adjust the selection 
criteria for the voxel selection algorithm until the points to be eliminated fall outside of 
the selection criteria for the displayed points that satisfy the voxel selection algorithm. 
Each tune the selection criteria is adjusted, the image must be interrupted. This iterative 
process is time consuming, and interferes with the visualization process of the user. 

Thus, there is a need in the art for a system and method for imaging 3D volume 
data sets that overcomes the deficiencies detailed above. Particxilarly, there is a need for 
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a system and method that re-draws images of large 3D volume data sets in response to 
user input at a rate sufficiently fast that the user perceives an instantaneous or real-time 
change in the image, wdthout perceptible delay or lag. There is a need for a system and 
method that allows a user to interactively change the displayed image in a contmuous 
manner, without interruption or perceptible delay or lag. Such a system and method 
would allow a user to more quickly and accurately interpret and ideadfy features inherent 
in 3D volume data sets. 
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Summary of the Invention 

The present invention is directed to a system and method for analyzing and 
imaging 3D volume data sets using a 3D sampling probe and other interactive tools. In 
one aspect of the invention, a method and computer program is provided which may be 
stored on a storage means for imaging a three-dimensional volume data set comprising 
a plurality of voxels where each voxel comprises a three-dimensional location and a 
15 dataword. The dataword may be representative of a physical phenomena, e.g., an 
amplitude signal related to a particular position within a geobody. The method may 
comprise steps such as creating at least one three-dimensional sampling probe, wherein 
the three-dimensional sampling probe is the same size or a subset of the three- 
dimensional volume data set. The three-dimensional sampling probe has a probe face 
20 plane and an opposing probe face plane. Other steps may include producing a pluraUty 
of control points in the probe face-plane such that the plurality of control points dejSne 
one or more lines on the probe fece plane, and extending a ribbon section 3&om the one 
or more Imes on the probe fece plane toward the opposing probe fece plane. An edge of 
the ribbon section is preferably formed by the one or more lines. Another step of the 
25 method may include selectively imaging datawords representative of the physical 
phenomena at three-dimensional locations which intersect the ribbon section and the 
three-dimensional sampling probe. 

The method may further comprise steps of editing the pluraUty of control points 
on the probe face plane to thereby redefine the one or more lines, and extending a 
30 correspondingly redefined ribbon section fiom the one or more lines on tiie probe fece 
plane toward the opposing probe face plane. The step of editing may fiirther comprise 
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functions such as deleting one or more of the plurality of control points, changing a 
location of one or more of the plurality of control points, and adding one or more control 
points to the piuralitj' of control points. 

In a prefeired embodiment, the ribbon section is perpendicular to the probe face 
plane and the ribbon section may extend from the probe face plane to the opposing probe 
face plane. The one or more lines forming the edge of the ribbon section may be edited 
through tlie plurality of control points to construct a plurality of open straight Imes or a 
closed line geometrical figure, if desired The ribbon section is preferably comprised of 
a plurality of planes. The ribbon section may or may not be parallel with respect to each 
of a plurality of side faces of the probe. 

In another embodiment of the present invention related to tracking a particular 
physical phenomena, such as a geological fault, the method may comprise the steps of 
positioning the probe face plane at a first position within the three-dimensional volume 
data set and forming a first set of control points on the probe face plane for tracking a 
physical phenomena described by the three-dimensional volume data set Another step 
may include interpolating between the first set of control points to define a first spline 
curve. Other steps may include moving the probe to a second position within the three- 
dimensional volume data set, forming a second set of control pomts on the probe face 
plane for tracking the physical phenomena and interpolating there between such that the 
second set of control pomts define a second splme curve. Another step may include 
inteipolating a three dimensional surface between the first spline curve and the second 
spline curve which is representative of the physical phenomena. 

The method fiirther permits displaying the interpolated surfece where the surface 
intersects the first set of control points and the second set of control points. It is an 
advantage of the present invention that the first splme curve, second splme curve and 
subsequent spline curves are curvilinear. 

Additional steps may include the reiterative process of moving the probe to a 
third position within the three-dimensional volume data set; forming a third set of control 
points on the probe face plane for tracking the physical phenomena, interpolating 
between the third set of control points to define a third spline curve, and interpolating 
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between the first spline curve, the second spUne, and the third spUne curve for fiirther 
defining the three dimensional surface representative of the physical phenomena. 

If desired, the method may fiirther comprise steps such as editing the 
respiesentive control points on the probe face plane at respective positions of the probe. 
Moreover, the method may include displaying a curvilinear connection ("v curves") 
between respective control points at respective positions of tlie probe. Another step may 
include displaying the spline curves and tiie v-curves on the three dimensional surface. 
The spline curves and the v-curves form a three dimensional grid also representative of 
tlie physical phenomena. The grid includes aplurality of mtersections between the spline 
curves and the v-curves. The method may fiirther comprise editing the current set of 
control points on the probe face plane, thereby reshaping the suifece and grid between 
the current spline curve and the prior spline curve. 

Preferablylhemefhod may also include steps such as selectmgoneofthepluraUty 
of mtersections to thereby reposition the probe fece plane to pass through the selected 
intersection. The metliod also comprises selecting one of the sets of control points to 
thereby reposition the probe face plane to pass through the selected set of control points. 

Stated another way, an embodiment of the metliod may comprise steps such as 
positioning the probe face plane at a plurality of positions within the three-dimensional 
volume data set, formmg a set of control points at each of the plurahty of probe fece 
plane positions such that each set of control points defines a related spline curve, 

repositioning thepiobefeceplanemdinterpolalingbetweentheplura%ofsplinec^^ 
to form a three dimensional surfece representative of the physical phenomena. 

Features and Advantages 

It is a feature of tiie present invention that a ribbon section tiirough a 3D sampling 
probe can be created, redrawn, edited, and moved quickly and convenientiy by creating 
aplurality of lines that are then projected through the 3D sampling probe. The lines may 
be drawn at angles offset fi-om flie coordinate system, such as an x, y. z or Cartesian 
coordinate system, of the 3D samplmg probe. 

It is another feature of the present invention tiiat structures in a 3D data volume 
set. such as for mstance geological stiuctures. can be quickly mapped by selecting points 



wo 02/37422 



PCT/USOO/29835 



of interest at a plurality of locations in the 3D sampling probe, which points may then be 
interpolated to produce a grid or surface related to the structure. The grid may be quickly 
edited and the probe may be moved to various points on the surface by selecting grid 
intersections. 

It is yet another feature of the present invention that, as.a user interactively moves 
a 3D sampling probe through a 3D volume data set, the image on the surfaces of the 3D 
sampling probe is re-drawn "on the fly" so that the user perceives the image changing in 
real-time with movement of the 3D sampling probe. Similarly, as a user interactively 
moves a 3D sampling probe through a 3D volume data set, the 3D sampling probe is 
volume rendered with varying degrees of transparency "on the fly" so that the user 
perceives the knage changing in real-time with movement of the 3D sampling probe. 

It is a further feature of the present mvention that a user can interactively change 
the shape or size of a 3D sampling probe so that the image on the surfaces of the 3D 
sampling probe is re-drawn "on tlie fly" so that the user perceives the image changmg in 
real-time with the change in shape or size of the 3D sampling probe. Similarly, a user 
can interactively change the shape or size of a 3D sampling probe so that the 3D 
sampling probe is volume rendered with varying degrees of transparency "on the fly" so 
that the user perceives the image changing in real-time with the change in shape or size 
of the 3D sampling probe. 

It is yet a further feature of the present invention that a user can interactively 
rotate a 3D sampling probe so that the image on the surfaces of the 3D sampling probe 
is le-drawn "on the fly" so that the user perceives the unage changing in real-time with 
the rotation of the 3D sampling probe. Similarly, a user can interactively rotate a 3D 
sampling probe so that the 3D sampling probe is volume rendered with varying degrees 
of transparency "on the fly" so that the user perceives the image changing in real-time 
with the rotation of the 3D sampling probe. 

It is yet a further feature of the present invention that an eraser 3D sampling probe 
can be created and manipulated by the user to directly delete from an unage selected 
points that fall within a certain datavalue range. 
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It is an advantage of the present invention that a user can manipulate a 3D 
sampling probe to interactively traverse a 3D volume data set to continuously foUow and 
im^e a feature. 

It is a further advantage of the present invention that a user can interactively 
change the displayed image in a continuous manner, without interruption or perceptible 
delay or lag. This aUows a user to more quicldy and accurately interpret and identify 
features inherent in 3D volume data sets. 

It is yet a further advantage of the present invention that the 3D sampling probes 
can be interactively re-shaped by the user to match the shape of geologic features, thereby 
enabling the user to better visualize and define the extent of geologic features. 

A still fijrther advantage of thepresent invention is that it can be used to visualize 
and interpret large volumes of 3D seismic data. The present invention can be used to 
quickly and accurately identify drilling sites. The present invention can advantageously 
be used to sharply reduce 3D seismic project cycle times, to boost production from 
existing wells, and to locate additional reserves. 

Brief Description of the Figures 

The present invention is described with reference to the accompanying drawings. 
In the drawings, like reference numbers indicate identical or fimctionally similar 
elements. AdditionaUy. the left-most dlgit(s) of a reference number identifies the 
drawing in which Hie reference number first appears. 

FIG. 1 shows one embodiment of a software or program structure for 

implementing the present invention; 
FIG. 2 shows a block diagram of one embodiment of the 3D sampling 

probe program of the present invention; 
FIG. 3 shows a curve illustrating opacity as a function of datavalue; 
FIG. 4 shows a flow diagram illustratmg one embodiment for 
implementing the present invention; 

FIG.5 showsaflowdiagramofoneembodimentforchangingadefeult 
probe; 
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FIG, 6 shows a flow diagram of one embodiment for creating 
additional probes; 

FIG. 7 shows a flow diagram of one embodiment for moving a probe; 
FIG. 8 shows a flow diagram of one embodiment for re-shaping a 
5 probe; 

FIG. 9 shows a flow diagram of one embodiment for rotating a probe 
in 3D space; 

FIG, 10 shows a flow diagram of one embodiment for rotating a probe 
while fixed in space; 

1 0 FIG. 11 shows a flow diagram of one embodiment for carrying out auto 

picking or seed selection; 
FIG. 12 shows one embodiment of a computer system suitable for use 

wifli the present invention; 
FIG. 13 shows an alternate embodiment of a computer system suitable 
15 for use with the present invention; 

FIG. 14 shows further detail of an exemplaiy computer system suitable 

for use with the present invention; 
FIG. 15 illustrates three opaque probes of the present invention, with 
two of the probes intersecting each other; 
20 HG. 16 illustrates three probes of the present invention, a data probe, a 

transparent cut probe, and a volume-rendered probe; 
FIG. 17 illustrates a ribbon section accordmg to the present mvention in 

tiie shape of a star; 
FIG. 18 shows a block diagram of a system for producing the ribbon 
25 section of FIG. 17; 

FIG. 19 illustrates a three dimensional grid and three dimensional 
surface representative of a physical phenomena described by a 
3-D volume data set according to the present mvention; and 
FIG. 20 shows a block diagram of a system for producing the grid and 
30 surface of FIG, 19. 
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Detailed DEscwraoN of the Preferred Embodiments 
Overview 

The present invention is directed to a system and method for analyzing and 
imaging three-dimensional ("3D") volume data sets using a 3D sampling probe. 3D 
volume data sets comprise "voxels" or volume elements. Each voxel is a sample or point 
within a volume. Each voxel can be expressed in the form (x, y, z, datavalue) where "x, 
y, z" identifies the 3D location of the point within the volume, and "datavalue" is tlie 
value of some measured or calculated attribute or physical parameter at the specified 
point within the volume. For example, a 3D volume data set suitable for use with the 
present invention is 3D seismic data. Each voxel in a 3D seismic data can be expressed 
as (x, y, z, ampUtude), with amplitude corresponding to the amplitude of reflected sound 
waves at the specified (x, y, z) location. 

Any form of information that can be e^qpressed m the voxel form (x, y, z, 
datavalue) is suitable for use with the present invention. In addition to seismic data, 
examples from the oil and gas industry include information from closely spaced well 
logs, gravity and magnetic fields, remote sensing data, and sidescan sonar image data. 
Other geologic or physical information could also include temperature, pressure, 
saturation, reflectivity, acoustical impedance, and velocity. 

Another application for the present invention is for mining. For example, the 
present invention can be used to visualize and interpret geologic and geophysical data to 
locate mining sites, to locate and track deposits to be mined, or to locate and track 
geologic features such as feults that would affect mining operations. The present 
invention also has application for clean up of toxic, hazardous, or other types of waste. 
For example, the present mvention can be used to visualize and inteipret data 
representing the geographic extent and distribution of the waste at a particular site. Such 
visualization and inteipretation is useful for prioritizing clean up at various sites, and for 
developing a clean-up plan for a particular site. 

The present invention can also be used with uif ormation outside of the oil and gas 
industry. For example, the present invention can be used for analyzing and imaging in 
the medical field, where the datavalue element of the voxel is obtained from a CAT 
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(computerized axial tomography) scanner or a magnetic resonance imaging (MRI) 
procedure. 

By way of explanation and example, the present invention will be described in 
detail below using 3D seismic data as the 3D volume data set. It is to be understood, 
however, that the present invention is not limited in any way to the use of 3D seismic 
data. 

The present invention is particularly useful as a visualization tool for interpreting 
3D seismic data. As used herein, the term "visualization" refers to the construction of a 
thiee-dimensional pictuie in the user's mind of physical or geologic features or physical 
parameters that are inherently present in 3D volume data sets. Such physical features or 
parameters are typically not apparent from conventional means of processing 3D data 
sets, such as scanning a series of cross-sections of tlie 3D volume data set, because of the 
mental reconstruction that needs to take place in order for a user to mentally "picture" the 
three-dimensional feature. Because of this mental reconstmction, it is diflScult to 
commxmicate and share among users the same 3D image. For example, the same 3D 
mental image of the terrain will not necessarily be present in the mind of every person 
that reads or analyzes a two-dimensional ("2D ") contour map of that terrain. Through the 
use of 3D computer graphics, users can visualize, and communicate and share, the same 
3D image of 3D volume data sets. By visualizing 3D seismic data, a team of geologists, 
geophysicists, and engineers can interpret the visualized data to make exploration and 
production decisions, such as drilling locations and well paths. 

To accompUsh the visualization fmiction,flie preset invention uses the c^ 
graphics techniques of texture moping and volume rendering. By "texture map" is 
meant wrapping or mapping a 2D picture or image onto a 2D or a 3D object. For 
example, a photograph of a person can be texture mapped onto a coffee cup. 

The term "volimie rendering" or "volume imaging" refers to drawing a three- 
dimensional object in a manner that conveys to a viewer the three-dimensional nature of 
the object, even though the viewer may be looking at a two-dimensional display or 
screen. Computer graphics technology makes use of coloring, lighting, and shading 
techniques to convey to the mind of the viewer what is high or low, behind or in front, 
light or dark, etc. The perspective or viewpoint can be changed so that the viewer can see 
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all sides of the 3D object. Volume rendering typically includes some type of 
transparency/opacity (opacity = 1 - transparency) control so that certain parts of the 3D 
object are more transparent, thereby allowing a viewer to "see through" outer surfaces of 
an object and view its internal structures. 

The present invention enables fast visualization and analysis of very large 3D 
volume data sets through the use of a "sampling probe", also referred to herein as a 
"probe" or "probe object". As explained in more detail below, the sampling probes of the 
present invention have numerous attributes, one of which is that tliey are typically created 
as a 3D sub- volume of tlie whole 3D volume data set to be visualized and analyzed. 

A number of sampling probes can be created, shaped, sized, and moved 
interactively by the user within the whole 3D volume data set The intersection of the 
sampling probe with the whole 3D volume data set is texture mapped onto the suifeces 
of tiie sampling probe, or volume rendered with varying degrees of transparency within 
the sampling probe. As used herein, the term "interactive" or "interactively" refers to 
changing or re-drawing an image m response to user input at a rate sufficiently fast that 
the user perceives an instantaneous or real-time change in the image, without perceptible 
delay or lag. In practice, a fiame rate of at least about 10 to 15 frames per second is 
sufficient to achieve interactive imaging as described herein. For example, as the user 
moves the sampling probe, such as by "cUddng and dragging" with a "mouse", the user 
pCTceives the texture on the surfeces of the sampling probe changing in "red-time" with 
movement of the sampling probe. As the sampling probe changes shape, size, orlocation, 
there is no user-perceivable delay or lag in imaging the texture, or, with varying degrees 
of transparency, volume-iendered attributes. In this manner, the user can interactively 
move the sampling probes through the whole 3D volume, and more easily and effectively 
visualize and interpret the features and physical parameters that are present within the 
geographic space represented by the whole 3D volume data set 

System Description 

The present invention may be implemented using hardware, software or a 
combination thereof, and may be implemented m a computer system or other processing 
system. One embodimrait of a software or program structure 100 for implementing the 
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present iiwention is shown in FIG. 1. At the base of program structure 100 is an 
operating system 1 02. Suitable operating systems 102 include, for example, the UNIX® 
operating system, or Windows NT® from Microsoft Corporation, or other operating 
systems as would be apparent to one of skill in the relevant art. 

Menu and windowing software 104 overlays operating system 102. Menu and 
windowing software 104 is used to provide various menus and windows to facilitate 
interaction with the user, and to obtain user input and instructions. Menu and windowing 
software 104 can include, for example, Microsoft Windows™, X Window System™ 
(registered trademark of Massachusetts Institute of Technology), and MOTIF™ 
(registered trademark of Open Software Foundation Inc.). As would be readily apparent 
to one of skill in the relevant art, other menu and windowing software could also be used. 

A basic graphics libraiy 1 06 overlays menu and windowing software 1 04. Basic 
graphics libraiy 106 is an application programming interface (API) for 3D computer 
graphics. The functions performed by basic graphics library 106 include, for example, 
geometric and raster primitives, RGB A or color index mode, display list or immediate 
mode, viewing and modeling transformations, lighting and shading, hidden surface 
removal, alpha blending (translucency), anti-aliasing, texture mapping, atmospheric 
effects (fog, smoke, haze), feedback and selection, stencil planes, and accumulation 
bxiffer. 

A particularly prefenred basic graphics library 106 is OpenGL®, available from 
Silicon Graphics, Inc. ("SGI"), Mountain View, California. The OpenGL® API is a 
multiplatform industry standard that is hardware, window, and operating system 
independent OpenGL® is designed to be callable from C, C++, FORTRAN, Ada and 
Javaprogrammmg languages. OpenGL® performs each of the fimctions listed above for 
basic graphics Ubrary 106. Some commands in OpenGL® specify geometric objects to 
be drawn, and others control how the objects are handled. All elements of the OpenGL® 
state, even the contents of the texture memory and the frame buffer, can be obtained by 
a client application using OpenGL®. OpenGL® and the client application may operate 
on the same or dififerent machines because OpenGL® is network transparent. OpenGL® 
is described in more detail in the OpenGL® Programming Guide (ISBN: 0-20 1 -63274-8) 
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and the OpenGL® Reference Manual aSBN: 0-201-63276-4), the entirety of both of 
which are incorporated herein by reference. 

Visual simulation graphics Ubrai3' 108 overlays basic graphics library 106. Visual 
simulation graphics library 108 is an API for creatmg real-time, multi-processed 3D 
visual simulation graphics applications. Visual simulation graphics library 1 08 provides 
functions that bundle together graphics library state control functions such as lighting, 
materials, texture, and transparency. These functions track state and the creation of 
display lists that can be rendered later. 

A particularly preferred visual simulation grapliics library 1 08 is IRIS Performer, 
available from SGI in Mountain View, California. IRIS Performer supports the 
OpenGL® graphics Ubrarydiscussedabove. IRIS Performer includes two main Ubiaries, 
libpf and libpr, and four associated libraries, libpfdu, libpfdb, Ubpfiii, and libpfutU. 

The basis of IRIS Performer is the performance rendering Ubrary Ubpr, a low- 
level Ubrary providing high speed rendering fiinctions based on GeoSets and graphics 
state control using GeoStates. GeoSets are collections of drawable geometry that group 
same-type graphics primitives (e.g., triangles or quads) into one data object. The GeoSet 
contains no geometry itself only pointers to data arrays and index arrays. Because all the 
primitives in a GeoSet are of the same type and have the same attributes, rendering of 
most databases is performed at maximum hardware speed. GeoStates provide graphics 
state definitions (e.g., texture or material) for GeoSets. 

Layeredabove Ubpris Ubp^ areal-time visual simulation aiviroiunent providing 
a high-performance multi-process database rendering system that optimizes use of 
multiprocessing hardware. The database utihty Ubrary, libpfdu, provides functions for 
defining both geometric and appearance attributes of 3D objects, shares state and 
materials, and generates triangle strips from independent polygonal inpuL The database 
library libpfdb uses the facilities of libpfdu, libp^ and libpr to import database files in 
a number of industry standard database formats. The libpfiii is a user interfece library 
that provides buUding blocks for writing manipulation components for user interfeces (C 

andC^H-progranmiing languages). FinaUy.theUbpfutil is the utiUty Ubrary thatprow^^ 
routines for implementing tasks such as Multichannel Option support and graphical user 
interfece (GUI) tools. 
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An application program that uses IRIS Performer and OpenGL® API typically 
carry out the following steps in preparing for real-time 3D visual simulation: 

1 . Initialize IRIS Performer; 

2. Specify number of graphics pipelines, choose the 
multiprocessing configuration, and specify hardware mode as 
needed; 

3 . Initialize chosen multiprocessing mode; 

4. Initialize frame rate and set frame-extend policy; 

5. Create, configure, and open windows as required; and 

6. Create and configure display channels as required. 

Once the application program has created a graphical rendering environment by 
carrying out steps 1 through 6 above, then the application program typically 
iterates through a main simulation loop once per frame. 

7. Compute dynamics, update model matrices, etc.; 

8. Delay until the next frame time; 

9. Perform latency critical viewpoint updates; 

10. Draw a frame. 

A 3D sampling probe program 110 of the present inventions overlays visual 
simulation graphics library 108. Program 110 interacts with, and uses the functions 
carried out by, each of visual simulation and graphics library 108, basic graphics library 
106, menu and windowing software 104, and operating system 102 in a manner known 
to one of skill in the relevant art. 

3D sampling probe program 110 of the present invention is preferably written in 
an object oriented programming language to allow the creation and use of objects and 
object fimctionality. A particularly preferred object oriented programming language is 
C++. In canyuig out the present invention, program 110 creates one or more probe 
"objects". As noted above, the probe objects created and used by program 110 are also 
referred to herein as sampling probes or probes. Program 110 manipulates the probe 
objects so that they have the following attributes. 
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A probe corresponds to asub-volume of alarger3D volume. Particularly, aprobe 
defines a sub-set that is less than the complete data set of voxels for a 3D volume data set. 
A probe could be configured to be equal to or coextensive with the complete data set of 
voxels for a 3D volume data set. but the functionality of the present invention is best 
carried out when the probe corresponds to a sub-volume and defines a sub-set that is less 
than thecompletedatasetofvoxelsforaSDvolume data set. For example, a 3D volume 
data set of seismic data can contain from about 500 MB (megabytes) to about 10 GB 
(gigab5^es) or more of data. A 2,500 square kilometer geographic space of typical 3D 
seismic data contains about 8 GB of data. A probe of the present invention for a 500 MB 
seismic data set would preferably contain about 10-20 MB of data. 

By using probes that are a sub-volume of the larger 3D volume, flie quantity of 
data that must be processed and le-drawn for each fiame of an image is dramatically 
reduced, thereby increasing the speed with which the image can be re-drawn. The 
volume of a three-dimensional cube is proportional to the third power or "cube" of the 
dimensions of the three-dimensional cube. Likewise, the quantity of data in a 3D volume 
data set is proportional to the third power or "cube" of its size. Tlierefore, the quantity 
of data in a sub-volume of a larger 3D volume will be proportional to the "cubed root" 
of the quantity of data in the larger 3D volume. As such, the quantity of data in a 
probe of thepresent invention will be proportional to the "cubed roof CV) of the quantity 
of data inlheSDvolumeofT^di it isasub-volume. By onlyhavingtoprocessliie sub- 
set of data that relates 10 the sub-volume of the probe, the present invention can re^w 
an image in response to user input at a rate sufficiently &stfhattiie user perceives an 
instantaneous or real-time change in the image, without perceptible delay or lag. 

The probes of the present invention can be interactively changed in shape and/or 
size, and interactively moved within the larger 3D volume. The outside geometry or 
surfaces of a probe can be interactively drawn opaque or texture mapped while the probe 
is being changed in shape and/or size or while the probe is being moved. The probe can 
be drawn or volume rendered with varying degrees of transparency while the probe is 

being changed in shape and/or size or moved, thereby revealmg the internal stmctures or 
features of the probe. 
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The 3D sampling probes of the present invention can have any shape, including 
rectangular shapes having one or more right angles and non-rectangular shapes having 
no right angles. The 3D sampling probes of the present invention can have orthogonal or 
perpendicular planes as outer surfeces (e.g., squares and rectangles), parallel planes as 
outer surfeces (e.g., parallelograms), or curved outer surfaces (e.g., spheres, ovals, or 
cylinders). The present invention is not limited to 3D sampling probes of any particular 
shape. The 3D sampling probes of the present invention can have arbitrary shapes, such 
as the shape of a geologic feature identified by a user. For example, as a user moves the 
3D samplmg probe through a 3D volume of seismic data, a geologic feature may be 
visualized and identified by the user. The 3D sampling probe can be interactively re- 
shaped by the user to match tiie shape of the geologic feature, thereby eaiabling the user 
to better visualize and define the exteait of that geologic feature. 

A probe can be used to cut into another probe, and the intersection of flie two 
probes can be imaged. A probe can be used to hi^ght data ui accordance with a seed 
selection algorithm, A probe can also be used to "erase" or delete data in accordance with 
a seed de-selection algorithm. These attributes will be explained in more detail below. 

FIG. 2 shows a block diagram of one embodiment of 3D samplmg probe program 
110. Program 110 mcludes a User Interfece Module (UIM) 210, a Graphics Processing 
Module (GPM) 220, and a Volume Sampling Module (VSM) 230. A 3D volume data 
set is illustrated as data volume 240, also referred to herein as a 3D volume. UIM 210 
and GPM 220 communicate via a bi-directional pathway 212. GPM 220 sends 
instructions and requests for data to VSM 230 via pathway 222. UIM 210 sends 
instructions and requests to VSM 230 via pathway 214. UIM 210 interacts with data 
volume 240 through pathway 216. 

Voxel data &om data volume 240 is transferred to VSM 230 via data pathway 
234. VSM230transfersdatatoGPM220viadatapathway232. Data volume 240 stores 
the 3D volume data set m a manner well known to one of skill in the relevant art. For 
example, the format for data volume 240 can consist of two parts, a volume header 
followed by the body of data that is as long as the size of flie data set. The volume header 
typically contains information, m a prescribed sequence, such as the file patii (location) 
of the data set, size, dimensions m the x, y, and z directions, annotations for the x, y, and 
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2 axes, annotations for the datavalue, etc. The body of data is a binary sequence of bytes, 
one or more bytes per data value, that can be ordered in the following manner. The fiist 
byte is the datavalue at volume location (x, y, z)=(0,0,0). The second byte is the 
datavalue at volume location (1,0,0), the thkd byte is the datavalue at volume location 
(2,0,0), etc. When the x dunension is exhausted, then the y dimension is incremented, 
and finally the z dimension is incremented. The present invention is not limited in any 
way to a particular data format for data volume 240. 

User Interface Module 210 handles the user interface to receive commands, 
instructions, and input data from the user. UIM 210 interfaces with the user through a 
variety of menus through which the user can select various options and settings, either 
through keyboard selection or through one or more user-manipulated input devices, such 
as a "mouse", or a 3D pointing device. UIM 210 receives user input as the user 
manipulates the input device to move, size, shape, etc. a 3D sampling probe. 

Hie primary functions carried out by UIM 210 will now be deswibed. UIM 210 
inputs from the user the identification of one or more 3D volume data sets (represented 
by data volume 240) to use for unaging and analysis. When a plurality of data volumes 
are used, the datavalue for each of the plurality of data volumes represents a different 
physical parameter or attribute for the same geographic space. By way of example, a 
plurality of data volumes could include a geology volume, a temperatme volume, and a 
water-saturation volume. The voxels in the geology volume can be ejqjressed inflie form 
(x, y, 2, seismic amplitude). The voxels in flie temperature volume can be expressed in 
fliefoim(x,y,2,'C). 

The voxds in Ae water-saturation volume can be expressed in the form (x, y, z, 
%saturation). The physical or geographic space defined by the voxels in each of these 
volumes is the same. However, for any specific spatial location (xo. yo, Zq), the seismic 
amplitude would be contained in the geology volume, the temperature in the temperature 
volume, and the water-saturation in the water-saturation volvime. 

UIM 210 inputs fi:om the vseac mformation to create one or more 3D samplmg 
probes. Such mformation mcludes size, shape, and initial location of the probe. Such 
mformation can also include imaging attributes such as color, Ughting, shading, and 
transparency (or opacity). By adjusting opacity as a fimction of datavalue, certain 
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portions of the data volume are more transparent, fliereby allowing a viewer to see 
through surfaces. An exemplary opacity curve 300 is shown in FIG. 3. Opacity curve 
300 illustrates opacity (1 -transparency) as a function of datavalue. As would be readily 
apparent to one skilled in the art, datavalues vrith greater opacity (less transparency) will 
mask the imaging or display of datavalues with lower opacity (more transparency). 
Conversely, datavalues will less opacity and greater transparency will permit the imaging 
or display of datavalues with greater opacity and lower transpai'ency. 

UIM 210 receives input from the user for sizing and shaping the 3D sampling 
probes. As described m more detail below, in a preferred embodiment of the present 
invention, the user changes the shape and/or size of aprobe by clicking onto "sizing tabs" 
on the probe, and making changes in the dimensions of the probe in one or more 
directions. UIM 210 receives input from the user to move the position or location of a 
3D sampling probe within the data volume. In a preferred embodiment, a user 
manipulates a mouse to "click" onto a surface of the probe to be moved, and then moves 
the mouse to move the probe throughout the geographic space defined by the data 
volume. 

UIM 210 receives input from the user to cany out "auto picking" processes. In 
an auto picking process, data points (voxels) are selected based upon a selection 
algorithm. In a prefened embodiment, the selection algorithm is based upon a seed point 
witiiin the 3D data volume. The selection algorithm then selects data points tiiat: (i) 
satisfy flie selection criteria or algorithm (e.g., have a datavalue within a specified filter 
range); and (ii) have a connectivity with or are connected to the seed point. Through 
UIM 210, the user is prompted to identify a seed point within the 3D volume, and to 
identify a filter range of datavalues used by the selection algorithm to "pick" the selected 
points. Preferably, the seed pomt is within one of the 3D sampling probes. 

UIM 210 also receives input from the user regardmg the content of the displayed 
image. For example, the user can preferably select the content of the displayed image. 
The content of the displayed im2^e could include only the 3D sampling probe, i.e., its 
intersection with the 3D volume. Additionally, the 3D sampling probe could be 
displayed either with or without a bounding box that defines the outer geometry of the 
probe. Alternatively, the displayed image could include the 3D sampling probe, as well 
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as die data that occupies the background xz, yz, and xy planes, and/or the data that 
occupies the 3D volume outside of the 3D sampling probe(s) being displayed. 

To cany out the foregoing functions, UIM 210 sends a request to Volume 
Sampling Module 230 to load or attach those 3D volume data sets identified by the user. 
UIM 210 communicates via pathway 212 with Graphics Processing Module 220 that 
carries out the display and imaging. 

The primary functions carried out by GPM 220 will now be described. GPM 220 
processes data for imaging of 3D sampling probes with the hghting, shading, 
transparency, and other attributes selected by the user. To do so, GPM 220 uses the 
functions avaUable through basic graphics Ubrary 106 and visual simulation graphics 
Ubrary 108 described above. The user can select (through UIM 210) to display only the 
one or more 3D sampling probes that have been created. Alternatively, the user can 
select to display one or more 3D sampling probes, as well as the 3D data volume outside 
of the probes, i.e. voxels within the 3D volume that do not intersect any of the 3D 
sampling probes that are being displayed. 3D sampling probes that are being displayed 
are referred to herein as "active probes". 

GPM 220 processes the re-shaping and move requests that are received by UIM 
210 from the user. GPM 220 draws the re-shaped 3D sampling probe in accordance with 
the user-selected attributes (color, lighting, shading, transparency, etc.). As the user 
20 inputs a change in shape for a 3D sampling probe, the image with selected attributes is 
re-drawn sufficiently fest to be perceived as real-time by the vs&l. Similarly, GPM 220 
draws the 3D sampling probe in the new position or location in accordance with the user- 
selected attributes (color, Kghting, shading, transparency, etc.). As the user moves the 
3D sampling probe through the 3D volume, the image of the 3D sampling probe with 
25 selected attributes is re-drawn sufficiently fast to be perceived as. real-time by the user. 

GPM 220 processes "auto picldng" requests that are received by UIM 210. GPM 
220 will image selected points widiin the 3D volume in accordance with the selection 
algorithm. Alternatively, GPM 220 wiU "erase" selected points within the 3D volume 
in accordance with the selection algorithm. 
^0 To cany out the foregoing functions, GPM 220 communicates via pathway 212 

with UIM 210 so that the information requested by the vsa is imaged or displayed with 
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the selected attributes. GPM 220 obtains the needed data from data volume 240 by 
sending a data request via pathway 222 to Volume Sampling Module (VSM) 230. 

The primary function of VSM 230 is to extract the appropriate data from data 
volume 240 at the request of GPM 220. VSM 230 receives requests for data from GPM 
220 via pathway 222, VSM 230 extracts the required data fiom data volume 240 and 
transfers die data to GPM 220 via data padiway 232 for processing and display. VSM 
230 also receives instructions from UIM 210 via pathway 214 to load or attach the 3D 
data volumes identified by the user. 

Turning now to FIG. 4, a flow diagram 400 illustrating one embodiment for 
implementing the present invention is shown. A start up or initialization process is 
shown in a step 402. In step 402, the user specifies the one or more data volumes (240) 
to be used. The specified 3D volume data sets are loaded from disk into main memory 
(a description of hardware suitable for canying out the present invention will be 
described in more detail below). A default 3D sampling probe is created, and drawn. 
The default 3D sampling probe is a sub-volume of the specified 3D volume(s) of 
arbitrary size and shape. The present invention is not limited to any particular size or 
shape for the default 3D sampling probe. 

By way of example of the present invention, the default 3D sampling probe can 
be a square (having equal dimensions in the x, y, and z directions). To draw the square 
default 3D sampling probe, the bounding geometry is first drawn with one edge of the 
bounding geometry located on the z axis. Data is then extracted from data volume 240 
by VSM 230 to draw the image of the intersection of the square default 3D sampling 
probe with the 3D volume (data volume 240). Particularly, data is extracted that 
corresponds to the intersection of the square default 3D sampling probe with the 3D 
volume in the xz, yz, and xy planes. This data is then sent by VSM 230 to GPM 220 so 
that it can be texture mapped onto the planes of the bounding box to provide an image 
of the square default 3D sampUng probe. 

In one embodiment of the present invention, the data that occupies the 
background xz, yz, and xy planes themselves, as well as the data that occupies the 3D 
volume outside of the default 3D sampling probe, are also imaged or displayed (in 
addition to tiie default 3D sampling probe) during start up step 402. Alternatively, start 
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up Step 402 can be carried out so that the data that occupies the background xz, yz, and 
xy planes, or the data that occupies (he 3D volume outside of the default 3D sampling 
probe, is not displayed or imaged. Preferably, the present invention is cairied out so that 
the user can selectively display, or not display, the data that occupies the background xz, 
yz, and xy planes, as well as the data that occupies the 3D volume outside of the active 
probes. 

In a step 404, UIM 210 is waiting to respond to user input or request. User input 
is received tlirough a user input device suitable for use with a computer, including but not 
limited to, a keyboard, mouse, joystick, trackball, rollerball, roller point, or other type of 
suitable pointing device, etc. Preferably, the user input device comprises a mouse or 
other similar device that enables the user to "cUck" on a particular displayed image, and 
"drag" that displayed image to another location. Such a user input device allows a user 
to move and re-shape displayed probes. Such a user input device also allows a user to 
activate drop-down menus, and to select the various options for the color, shading, 
Ughting, and transparency attributes. A keyboard can also be used for entering 
information relating to the selected attributes. 

Reference numeral 406 refers generally to a piuraUty of functions that can be 
carried out by the present invention. These fimctions can be cairied out individually or 
smiultaneously,dependmguponmputfiiom the user. For example, aprobecanbemoved 
(fimction430)androtated(function450)simultaneously. While the fimctions identified 
by reference numeral 406 are beiig cairied out, the image of the 3D samphng probes is 
being re-drawn suffidentiy fest to be perceived as real-time by the user. Each of the 
functions identified by referttice numeral 406 will now be described. 

If a user wants to change the default probe, then function 410 is carried out. The 
steps for carrying out fimction 410 are shown m FIG. 5 by way of flow diagram 
connector 5 A. In a step 502, the changes to the defeult probe are input by UIM 210 from 
the user. For example, the changes to the defeult probe can be to the shape or size, the 
location, or the attributes such as color, shadmg, lighting, and transparency. 

In a step S04, UIM 210 sends a request to GPM 220 to draw the changed defeult 
probe. In a step 506, GPM 220 requests data for the changed defeult probe ftom VSM 
230. In making this request, GPM 220 would invoke function 430 if it was necessary to 



-22- 



wo 02/37422 



PCT/USOO/29835 



move the default probe, function 440 to re-shape the default probe, and functions 450 or 
460 to rotate the default probe. The foregoing functions will be described in more detail 
below. 

The data that will be extracted from data volume 240 by VSM 230 in response 
to the request made by GPM 220 in step 506 will depend upon attributes that have been 
selected by the user. If the opacity settings selected by the user are such that all 
datavalues are opaque, then the data extracted by VSM 230 will be limited to the surfaces 
of the changed default probe. Because of tlie selected opacity, it will not be possible for 
the user to see inside the changed default probe, so only the data corresponding to the 
surfaces or outside of the changed default probe will be extracted by VSM 230. In a step 
508, GPM 220 processes the data extracted by VSM 230 for tiie surfaces of tiie changed 
default probe, and draws the changed default probe by texture mapping onto the surfaces 
in accordance with the attributes selected by the user. By extracting only the data that 
can be seen by the user, the image of the changed default probe can be dmwn more 
quickly because less data needs to be processed, i.e., the data corresponding to the 
"inside" of the changed default probe is not processed. 

Alternatively, if the opacity settings selected by the user are such that some of tlie 
datavalues are opaque and some of the datavalues are transparent, then the data extracted 
by VSM 230 will include the data corresponding to the entire volume of the changed 
default probe. Because of the selected opacity and transparency, it will be possible for the 
user to see inside the changed defiult probe, so data cotresponding to the entire volume 
of the changed default probe will be extracted by VSM 230, In such a situation, GPM 
220 processes the data extracted by VSM 230 in step 508, and dmws the changed default 
probe by volume rendering in accordance with the attributes selected by the user. 

If a user wants to create additional probes, then function 420 is carried out. The 
present invention is not limited to any particular number of active probes. The steps for 
canying out function 420 are shown in FIG. 6 by way of flow diagram connector 6A. In 
a step 602, the shape, size, location, attributes, etc. for the additional probes are input by 
UIM 210 from the user. In a step 604, UIM 210 sends a request to GPM 220 to draw the 
additional probes. 
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In a Step 606, GPM 220 requests data for ttie additional probes from VSM 230, 
In a manner similar to that described above for changing the default probe, the data that 
is extracted from 3D or data volume 240 by VSM 230 will depend upon the opacity 
selected by the user for the additional probes. If the opacity settings selected by the user 
are such that all datavalues for the additional probes are opaque, then the data extracted 
by VSM 230 will be limited to the surfaces of the additional probes. Alternatively, if the 
opacity settings selected by tlie user for the additional probes are such that some of the 
datavalues are opaque and some of tiae datavalues are transparent, then tlie data extracted 
by VSM 230 will include the data corresponding to the entire volumes of the additional 
probes. In this manner, the additional probes can be drawn more quickly by minimizing 
the quantity of data Aat must be processed. 

In a step 608, GPM 220 processes the data extracted by VSM 230 for the 
additional probes, and draws the additional probes in accordance with the attributes 
selected by the user, either by texture mapping onto the surfaces of the additional probes, 
or by volume rendering the entire volume of the additional probes. 

If a user wants to move a probe, dien function 430 is carried out. The steps for 
carrying out fimction 430 are shown in FIG. 7 by way of flow diagram connector 7A. In 
a step 702, the new location for the probe is input by UIM 210 from the user. In a 
preferred embodiment of the present invention, the user inputs the new location of the 
probe by clicking a mouse or other type of suitable user input device to snap a pointer 
onto a surfece of the probe to be moved The user changes the location of the probe by 
movir^ the mouse or other suitable user input device in any direction, thereby dragging 
the probe along a trajectory. 

In a step 704, UIM 210 sends a move request to GPM 220 to draw the probe at 
the new location. GPM 220 requests data for the new location of the probe from VSM 
230. In a manner similar to that described above, the data that is extracted from data 
volume 240 by VSM 230 will depend upon the opacity selected by the user for the probe 
being moved. If the opacity settings selected by the user are such that all datavalues for 
the probe being moved are opaque, then the data ex^tracted by VSM 230 will be limited 
to the surfaces of the probe being moved. Alternatively, if the opacity settings selected 
by the user for the probe being moved are such that some of the datavalues are opaque 
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and some of the datavalues are transparent, then the data extracted by VSM 230 will 
include the data corresponding to the entire volume of the probe being moved. In this 
manner, the probe can be drawn at its new location more quickly by minimizing the 
quantity of data that must be processed. 

In a step 708, GPM 220 processes the data extracted by VSM 230 for the probe 
being moved, and draws the probe at its new location in accordance with the attributes 
selected by the user, either by texture mapping onto the surfaces of the probe being 
moved, or by volume rendering the entire volume of the probe being moved. 

As the user moves the probe, for each new location of the probe, steps 702 
through 708 are repeated at a rate sufficiently fast that the user perceives the image of the 
probe, with texture mapping or volume rendering as appropriate, changing in "real-time" 
with movement of the probe. The image is being re-drawn at a frame rate sufficiently 
fast to be perceived as real-time by the user. 

If a user wants to re-shape a probe, then function 440 is carried out As used 
herein, the term "re-shape" refers to any change in dimension of a 3D samplmg probe in 
any direction. The shape of a 3D sampling probe can be changed, or re-shaped, for 
example, by changing the size in one or more directions, such as by changing a square 
probe into a rectangular probe by increasing the size of the probe in the x direction, and 
decreasing the size of the probe in the y direction. As another example, the shape of a 
3D sampling probe can be changed by changmg the siu^e from spherical to rectangular. 
As yet anotiier example, a square 3D sampling probe (equal dimensions in the x, y, and 
2 duections) can be re-shaped in accordance with the present invention to be a larger or 
smaller square-shaped probe by changing the size equally in each of the x, y, and z 
directions. The re-shaped probe also has a square shape, but as a larger or smaller square. 

The steps for carrymg out function 440 are shown in FIG. 8 by way of flow 
diagram connector 8A. In a step 802, the new shape and/or size for the probe is input by 
UIM 210 from the user. In a preferred embodiment of the present mvention, the user 
inputs the new shape of a probe by clicking a mouse or other type of suitable user input 
device to snap a pomter onto a "sizmg tab" of the probe to be re-shaped. As used herein, 
a "sizing tab" refers to a designated area on a surface of the probe. Such a designated 
area is preferably displayed in a color that is different from the colors bemg used to 
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display the features or physical parameters of the 3D volume data set. When the pointer 
is snapped to the sizing tab, manipulation of the mouse or user, input device changes the 
dimensions or proportions of flie surface on w*ich the sizing tab is located. When the 
desired size or shape is reached, the user again cUcks the mouse or user input device to 
release the pointer from the sizing tab. Sizing tabs are illustrated in FIGS. 1 5 and 1 6. The 
sizing tabs are the small dark squares that appear on the surfaces of the probes, along the 
bounding geometry of the probes. Hie location of the sizing tabs is not limited to tlie 
bounding geometry of the probes. The user changes the shape of the probe by clicking 
the mouse or other suitable user input device onto a sizing tab. moving the mouse until 
the surface being changed has the desired shape, and then releasing the mouse from the 
sizing tab. This process can be repeated, if necessary, using other sizing tabs on the 
probe until the probe is re-shaped to the desired shape. 

It would be readily apparent to one of skill in the relevant art how to implement 
such a sizing tab for re-shaping the probes of the present invention. It is to be 
understood, however, that the present mvention is not limited to the use of sizing tabs for 
re-shaping probes, and other suitable methods can be used. For example, the user could 
select from a number of pre-set sliapes (e.g., squares, rectangle, cylinders, spheres) by 
activating a drop-down menu, or by scrolUng through the shapes by repeatedly clicking 
a mouse. 

In a step 804, UIM 210 sends a re-shape request to GPM 220 to draw the re- 
shaped probe. In a step 806, it is detennined whedier more data is needed to draw the re- 
sh^)ed probe. For exan^le, if the re-shaped probe is of a diape and size fliat "fits inside" 
the existing probe, tiien no more data is needed, and processing continues at a step 810. 
Alternatively, if the re-shaped probe is of a shape and size that falls at least partiaUy 
outside of the existing probe, then, in a step 808, GPM 220 requests the data needed for 
the re-shaped probe from VSM 230. hi a manner sunUar to that described above, the data 
tliat is extracted from 3D or data volume 240 by VSM 230 will depend upon the opacity 
selected by the user for the probe being re-shaped. If die opacity settings selected by flie 
user are such that all datavalues for flie probe being re-shaped are opaque, then flie data 
extracted by VSM 230 will be limited to the surfaces of the probe bemg re-shaped. 
Alternatively, if the opacity settings selected by the user for the probe being re-shaped 
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are such that some of the datavalues are opaque and some of the datavalues are 
transparent, then the data extracted by VSM 230 will mclude the data corresponding to 
the entire volume of the probe being re-shaped. In this manner, the probe can be drawai 
with its new shape more quickly by minimizing the quantity of data that must be 
processed. 

In step 810, GPM 220 processes the data extracted by VSM 230 for the probe 
being re-shaped, and draws the probe with its new shape in accordance with the attributes 
selected by the user, either by texture mapping onto the surfaces of tlie probe being re- 
shaped, or by volume rendering the entire volume of the probe being re-shaped. 

As the user changes the shape of the probe, steps 802 through 810 are repeated 
at a rate suflSciently fast that the user perceives the image of the probe, with texture 
mapping or volume rendering as appropriate, changing in *'real-time" with the changing 
shape of the probe. The image is being re-drawn at a frame rate sufficiently fast to be 
perceived as real-time by the user. 

If a user wants to rotate a probe in 3D space, then function 450 is carried out. In 
function 450, tlie 3D orientation, wliich is the same for both the 3D volume and the 
probe, is changed, thereby rotating the 3D volume and the probe in space. The steps for 
carrying out fiinction 450 are shown in FIG. 9 by way of flow diagram connector 9A. In 
a step 902, the new 3D orientation for the 3D volume and the probe is input by UIM 210 
ftom the user. In a preferred embodiment of the present invention, the user inputs the 
new orientation by clicking a mouse or other ^rpe of suitable user input device to snap 
apointer onto an axis of the probe to be rotated. Manipulation of flie mouse or user input 
device changes the orientation of that axis. When the desired orientation is reached, the 
user again clicks the mouse or xiser input device to release tiie pomter from tixe axis. It 
would be readily apparent to one of skill in the relevant art how to implement such a 
change in orientation. It is to be understood, however, that the present invention is not 
limited to changing the orientation in this manner. For example, the user could select 
from a number of pre-set rotations (e.g., rotate 90' to the left or right; rotate 45' to the left 
or right, etc.) by activating a drop-down menu, or by scrolling through the rotations by 
repeatedly clicking a mouse. 
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In a st^ 904, UIM 210 sends a request to rotate in 3D space to GPM 220 to draw 
the rotated probe. Li a step 906, GPM 220 requests data for the rotated probe from VSM 
230. In a nianner similar to that described above, the data that is extracted from 3D or 
data volume 240 by VSM 230 will depend upon the opacity selected by the user for the 
probe being rotated. If the opacity settings selected by the user are such that all 
datavalues for the probe being rotated are opaque, then the data extracted by VSM 230 
will be limited to the surfaces of the probe being rotated. Alternatively, if the opacity 
settings selected by the user for the probe being rotated are such that some of the 
datavalues are opaque and some of the datavalues are transparent, then the data extracted 
by VSM 230 will include the data corresponding to the entire volume of the probe being 
rotated. In this manner, the probe can be drawn with its new orientation more quickly by 
m i nim izin g the quantity of data fliat must be processed. 

In step 908, GPM 220 processes the data extracted by VSM 230 for the probe 
being rotated, and 6mvrs the probe with its new orientation in accordance wifli the 
attributes selected by the user, either by texture mapping onto the surfaces of the probe 
being rotated, or by volume rendering the entire volume of the probe being rotated. 

As the user rotates the probe in 3D space, steps 902 through 908 are repeated at 
a rate sufficiently fest that the user perceives the image of flie probe, with texture 
mapping or volume rendering as appropriate, changing in "real-time" with the changing 
orientation of the probe. The image is being re-drawn at a frame rate sufficiently fast to 
be perceived as real-time by the itser. 

If a user wants to rotate a probe while it is fixed in 3D space, then function 460 
is carried out In function 460, the 3D orientation of the probe is rotated independentiy 
of the 3D orientation of the 3D volume, thereby rotating the probe while it is fixed in the 
3D space defined by die orientation of the 3D volume. In this manner, the background 
planes for an active probe can be displayed in a fixed orientation, and flie active probe 
can be rotated witiiin the back^ound planes. 

The steps for carrymg out function 460 are shown in FIG. 10 by way of flow 
diagram connector 1 OA. In a step 1002, the new 3D orientation for the probe is input by 
UIM 210 from the user. In a preferred embodiment of the present mvention, the user 
selects the option to rotate wtdlc fixed in space, for example, from a "drop-down" menu. 
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Tlie user then inputs the new orientation for the probe by clicking a mouse or other type 
of suitable user input device to snap a pointer onto an axis of the probe to be rotated. 
Manipulation of the mouse or user input device changes the orientation of that axis. 
When the desired orientation is reached, the user again clicks the mouse or user input 
device to release the pointer from the axis. It would be readily apparent to one of skill 
in the relevant art how to implement such a change in orientation. It is to be understood, 
however, that the present invention is not lunited to changing tlie orientation in this 
mamier. For example, the user could select from a number of pre-set rotations (e.g., 
rotate 90' to the left or right; rotate 45^ to the left or right, etc.) by activating a drop-down 
menu, or by scrolling through the rotations by repeatedly cUcking a mouse. 

In a step 1004, UIM 210 sends a request to rotate while fixed in space to GPM 
220 to draw the rotated probe. In a step 1006, GPM 220 requests data for the rotated 
probe from VSM 230. In a manner similar to that described above, tlie data that is 
extracted from 3D or data volume 240 by VSM 230 will depend upon the opacity 
selected by the user for the probe being rotated. If the opacity settings selected by the 
user are such that all datavalues for tlie probe being rotated are opaque, then the data 
extracted by VSM 230 will be limited to the surfaces of the probe being rotated. 
Alternatively, if the opacity settings selected by the user for the probe being rotated are 
such that some of the datavalues are opaque and some of flie datavalues are transparent, 
then the data extracted by VSM 230 will include the data coirespondmg to the entire 
volume of the probe being rotated; In fins manner, the probe can be drawn with its new 
orientation more quickly by minimizing the quantity of data fbat must be processed. 

In step 1008, GPM 220 processes the data extracted by VSM 230 for the probe 
being rotated, and draws the probe with its new orientation in accordance with the 
attributes selected by the user, either by texture mapping onto the surfaces of the probe 
being rotated, or by volume rendering the entire volume of the probe being rotated. 

As the user rotates the probe while it is fixed in space, steps 1002 through 1008 
are repeated at a rate suflSciently fast that the user perceives the image of the probe, with 
texture mappijD^ or volume rendering as appropriate, changing in "real-time" with the 
changing orientation of the probe. The image is being re-drawn at a fi:ame rate 
sufiBciently fast to be perceived as real-time by the user. 
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If a user wants to cany out an "auto picking" process, then function 470 is cairied 
out. The steps for carrying out function 470 are shown in FIG. 11 by way of flow diagram 
connector 1 1 A. In a step 1102, a seed point within the data set of the 3D volume, and a 
selection criteria based on datavalues, are input by UIM 210 &om the user. Preferably 
tlie seed point is within the data set of voxels that defines a probe. As described below, 
such a probe is referred to herein as a seed 3D sampling probe or aii eraser 3D sampling 
probe. However, the seed point can be witliin the data set of voxels that defines the 3D 
volume, outside of an active probe. In a preferred embodiment of the present invention, 
the user selects the option to execute an auto picking process, for example, from a "drop- 
down" menu. The user then selects the seed point by clicking a mouse or other type of 
smtable user input device to snap a pointer onto the desired seed point. The selection 
criteria can be input, for example, by graphically selecting a range, or by keying in 
specific numerical values. It would be readily apparent to one of skUl in the relevant art 
how to input firom the user a seed point and filter range of datavalues. 

In a step 1104, UIM 210 sends an auto picking request to GPM 220 to draw the 
rotated probe. In a step 1106, GPM 220 requests selected points to image from VSM 230. 
The selected points are those tliat are connected to the seed point, and that have a 
datavalue within the selection criteria. 

In step 1108, GPM 220 processes the data extracted by VSM 230 to draw the 
selected points. The selected points are preferably highfighted by being drawn in a color 
diffaent from those used to depict the features or physical parameters of the 3D volume 
data set Alternatively, step 1108 can be earned out to "erase" or delete from the image 
the selected points. 

In a similar manner, auto picking function 470 can be used to "erase" or de-select 
points. For example, an eraser 3D sampling probe is defined, such as by invoking 
fimction 420 to create an additional probe. A "de-selection" criteria based on datavalues 
is defined. Points previously selected by an auto picking operation that satisfy the de- 
selection criteria are identified as candidates for de-selection. As the eraser 3D sampling 
probe moves through the 3D volume, the de-selected points are deleted from the image, 
and the image is re-drawn sufBciently fast to be perceived as real-time by the usex. 
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Once auto picking function 470 is initiated by the user, it can be carried out 
simultaneously with, for example, move fimction 430. In this manner, as the user moves 
the probe, steps 1102 through 1108 (and steps 702 through 708) are repeated at a rate 
sufficiently fast that the user perceives the image of the probe, with the selected points, 
changuxg in "real-tune" with the changing location of the probe. As the probe is moved, 
the selected points can be highlighted by being drawn in a suitable color, thereby having 
the auto-piclcing 3D sampling probe function as a "highlighter" as it moves through the 
3D volume. Alternatively, as the probe is moved, points previously selected by an auto 
picking operation can be "erased" or deleted from the image, thereby having the probe 
function as an "eraser" or eraser 3D sampling probe as it moves through the 3D volume. 
In either embodiment, the image is being re-drawn at a frame rate sufBciently fast to be 
perceived as real-time by the user. 

If a user wants to create a "ribbon section," then fimction 480 is performed. Tlie 
steps necessary for performing function 480 are described further below in reference to 
FIG. 17 and by way of block diagram connector 18A in FIG. 18. 

If a user wants to create a "3D surface' representative of a physical phenomena 
found widiin a 3D volume data set, then fimction 490 is performed. The steps necessary 
for performing function 490 are described further below in reference to FIG. 19 and by 
way of block diagram connector 20A in FIG. 20. 

In any event where a user desires to cany out one or more of the functions 
described above such as more probe (430), re-shape piobe (440), create a ribbon section 
(480) and create a 3D surface (490), each function can be performed independent of; or 
in connection with, one or more of the other functions. 

With reference now to FIG. 12, one embodiment of a computer system suitable 
for use with the present invention is shown. A graphics supercomputer 12 10 contains one 
or more central processing units (CPU) or processors 1212. Supercomputer 1210 
contains a random access memory (RAM) 1214 that can be accessed by processors 1212. 
Supercomputer 1210 also contains one or more graphics modules 1216 that also access 
RAM 1214. Gr^hics modules 1216 execute the functions carried out by Graphics 
Processing Module 220, using hardware (such as specialized graphics processors) or a 
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combination of hardware and software. A user input device 1218 allows auser to control 
and input information to graphics supercomputer 1210. 

A particularly preferred graphics supercomputer is an Onyx2 Infinite Reality 
system, availablefi-om SUicon Graphics, Inc., Mountain View, CA, configured with eight 
processors, three graphics pipeUnes, 16 GB of main memory, and 250 GB of disk 
memory. Such a graphics supercomputer has a scalable, high-bandwidth, low-latency 
■architectui-e to provide high speed reuderiug on multiple graphics pipelines. Graphics 

supercomputers fromotliervendors, such as Hewlett-Packard Company ofPalo Alto, CA 
or Sun Microsystems of Mountain View, CA could also be used. 

The graphics data forming the image to be displayed is sent from graphics 
supercomputer 1210 to a multiple-screen display system 1220 for projection onto a 
screen 1230. In the embodiment shown in FIG. 12, three projectors are used. From the 
perspective of a user viewing the image on screen 1230, the three projectors include a left 
projector 1240, a center projector 1250, and a right projector 1260. Although three 
projectors are shown, the present invention is not limited to tlie use of any particular 
number of projectors. 

Projector 1240 has a projection field on screen 1230, shown generally at 1242, 
between a point 1241 and a point 1243. Projector 1250 has a projection field on screen 
1230, shown generally at 1252, between a point 1251 and a point 1253. Projector 1260 
has a projection field on screen 1230. shown generally at 1262, between a point 1261 and 
a point 1263. Projection fields 1242 and 1252 have an overly region 1244, between 
points 1251 and 1243. Similariy, projection fields 1262 and 1252 have an overlap region 
1264, between points 1261 and 1253. The hnage to be displayed is divided into three 
(left, center, and right) over-lapping sub-images. By simultaneously projecting the thi«e 
over-lapping sub-images, thefield-of-view to theuser is increased over that avaUable, for 
example, on a monitor or through the use of just one projector. As an example, use of the 
three over-lapping sub-images shown in FIG. 12 incieases the field-of-view to 
approximately 160*. Overlap regions 1244 and 1264 are each approxunately 5.3'. 
Multiple-screen display system 1220 accounts for overlap regions 1244 and 1264 in a 
weU-known manner to edge-blend the images of the three projectors to form one 
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seamless image on screen 1230. Suitable display and projector systems are available 
from SEOS, London, England, such as the Barco projector units. 

FIG. 13 shows an alternate embodiment of a computer system suitable for use 
with the present uivention. In the embodiment shown in FIG. 13, graphics 
supercomputer 1210 is configured with multiple processors 1212, RAM 1214, and two 
graphics modules 1216. Graphics workstations suitable for use in the embodiment shown 
in FIG. 1 3 are available from Silicon Graphics, Inc. or Sun Microsystems. Each graphics 
module 1216 is connected to a monitor 1320 for display. Monitor 1320 should preferably 
be a color graphics monitor suitable for display of grapliics such as that shown in FIGS. 
15 and 16. Preferably, one of monitors 1320 displays the image of the 3D sampling 
probes, and the otiier of monitors 1320 displays the various menus used to operate 3D 
samplmg probe program 1 10. FIG. 13 also shows a keyboard 1330 and a mouse 1332 
that fimction as user input devices. 

A computer system capable of carrying out the fiinctionality described herein is 
shown in more detail in FIG. 14, Computer system 1402 includes one or more 
processors, such as processor 1404. Processor 1404 is connected to a commmiication bus 
1406. Various software embodiments are described in terms of this exemplary computer 
system. After reading this description, it will become apparent to a person skilled in the 
relevant art how to implement the invention using other computer systems and/or 
computer architectures. 

Computer system 1402 also includes a main memory 1408, preferably random access 
memory (RAM), and can also include a secondary memory 1410. Secondary memory 
1410 can include, for example, a hard disk drive 1412 and/or a removable storage drive 
1414, representing a floppy disk drive, a magnetic tape drive, an optical disk drive, etc. 
Removable stor^e drive 1414 reads from and/or writes to a removable storage unit 14 18 
in a well known manner. Removable storage unit 1418, represents a floppy disk, 
magnetic tape, optical disk, etc. which is read by and written to by removable storage 
drive 1414. As will be appreciated, removable storage unit 1418 includes a computer 
usable storage medium havmg stored therein computer software and/or data. 

In alternative embodiments, secondary memory 1410 may include other similar 
means for allowing computer programs or other instmctions to be loaded into computer 
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system 1402. Such means can include, for example, a removable storage unit 1422 and 
an interfece 1420. Examples of such can mclude a program cartridge and cartridge 
interface (such as that found in video game devices), a removable memory chip (such as 
an EPROM, or PROM) and associated socket, and other removable storage units 1422 
and interfaces 1420 which allow software and data to be transferred from removable 
storage unit 1422 to computer system 1402. 

Computer system 1402 can also include a communications interface 1424. 
Communications interface 1424 allows software and data to be transferred between 
computer system 1402 and external devices. Examples of communications interface 
1424 can include a modem, a network interface (such as an Ethernet card), a 
communications port, a PCMCIA slot and card, etc. Software and data transfened via 
communications interface 1424 are in the form of signals 1426 that can be electronic, 
electromagnetic, optical or oflier signals capable of being received by communications 
interface 1424. Signals 1426 are provided to communications interface via a channel 
1428. Channel 1428 carries signals 1426 and can be implemented using wire or cable, 
fiber optics, a phone line, a ceUular phone link, an RF link and other communications 
channels. 

In this document, flie terms "computer program medium" and "computer usable 
medium" are used to generally refer to media such as removable storage device 1418, a 
hard disk installed in hard disk drive 1412, and signals 1426. These computer program 
products are means for providing software to computer system 1402. 

Computer programs (also called computer control logic) are stored in main 
memory 1 408 and/or secondary memory 1410. Computer programs can also be received 
via communications interfece 1424. Such computer programs, when execute4 enable 
computer system 1402 to perform the features of the present invention as discussed 
herein. In particular, the computer programs, when execute4 enable processor 1404 to 
perform the features of the present invention. Accordmgly, such computer programs 
represent controlIa:s of computer system 1402. 

In an embodiment where flie invention is implemented using software, the 
software may be stored in a computer program product and loaded into computer system 
1402 usuig removable storage drive 1414, hard drive 1412 or communications interface 
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1424. The control logic (software), when executed by processor 1404, caiises processor 
1404 to perform the functions of the invention as described herein. 

In another embodiment, the invention is implemented primarily in hardware 
using, for example, hardware components such as application specific integrated circuits 
(ASICs). Implementation of such a hardware state machine so as to perform the 
functions described herein will be apparent to persons skilled in the relevant art(s). 

In yet another embodiment, tlie invention is implemented using a combination of 
both hardware and software. 

System Operation and Results 

The operation and results of the present invention will now be described, usir^ 
a data volume 240 that contains seismic data (datavalues representing seismic 
amplitudes). The user specifies the particular seismic data volume to be used, which is 
loaded from disk into main memory. A default 3D sampling probe is drawn. The user 
specifies the colors to be used for tbe seismic amplitudes. The degree of transparency can 
also be selected. The tliree probes shown in FIG. 15 are all opaque, with the intersection 
of the probes and the seismic data volume texture mapped onto the surfaces of the probes. 
One of the probes is displayed with the bounding geometry showr^ the other two probes 
are displayed without the bounding geometry. 

FIG. 1 5 shows three active probes. The user has selected not to display the data 
contained in the background planes and in flie remainder of flie seismic data volume 
outside of the active probes. Two of the probes shown in FIG. 15 mteisect each other, 
and the intersection of the two probes is displayed. In this maimer, the user can more 
readily visualize and mterpret geologic features inherent in the seismic data volume. For 
example, a geologic feature, represented by a dark band between tWo light bands, extends 
across the face of the larger intersecting probe and "turns the comer" to extend onto the 
face of the smaller intersecting probe perpendicular to it. The ability to move and 
intersect the probes with each other, and throughout the seismic data volimie, enables a 
user to better interpret and track the extent of such a geologic feature. 

FIG. 1 6 illustrates how one probe can be used to "cut" another probe to create a 
"hole" in a probe. As with FIG. 1 5, the user has selected not to display the data contained 
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in the background planes and in fee remainder of the seismic data volume outside of tfie 
active probes. FIG. 16 illustrates that Ae opacity settings can be individually selected by 
the user for each active probe. One of the probes shown in FIG. 1 6 is opaque, so that it 
is not possible to see through the surfeces of this probe. In order to see internal to this 
outer probe, it must be cut away by anoflier probe. The outer opaque probe is referred 
to as a "data probe". A second completely transparent "cut probe" has been used to cut 
out a 3D sub-section of the data probe. Because the cut probe is completely transparent, 
it is not visible in FIG. 16. However, the fact that the completely transparent cut probe 
is present is evidenced by the fact that the opaque internal surfaces of the data piobe are 
visible. The image of the intersection of the data probe and the cut probe is the 
intersecting sur&ce internal to the data probe. 

A third active probe is shown in FIG. 1 6. The third probe is displayed with the 
bounding geometry shown. The third probe is volume rendered with varying degrees of 
transparency so that the user can see through the outer surfaces of the probe and view 
geologic features within the third probe. As shown in FIG. 1 6, the third probe is volume 
rendered partially within the 3D sub-section of the data probe that has been cut away by 
the cut probe. 

The third volume-r«idered probe shown in FIG. 1 6 also contains selected points 
that have been selected through a seed picking process (function 470). TTie selected 
points have been iniaged in a manner to highlight them for the user . The selected points 
are shown in FIG. 16asconnectingpoints. The seed point is illustrated in FIG. 16bythe 
darker sphere. 

FIG. 17 illustrates one embodiment of a ribbon section 1710 in accord with the 
present invention. Like the probe, ribbon sections are a 3-D volume visualization metitod 
for displaying data along user-defined traverses through a 3-D volume data set within a 
probe. The traverses cut through the 3-D volume data set somewhat like a cookie cutter 
through dough and therefore are referred to herein as cookie planes such as cookie planes 
1712 and 1714. Ribbon sections may display the 3-D data in orientations that do not 
necessarily conform to the orientation of the ordinate axis of the 3-D volume data set 
and/or probe. Ribbon section 1710 is produced by the user in accord with the process 
describedherembelow. The traverses or cookie planes are defined by digitizmg control 
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points such as control points 1716, 1718, and 1720, which may be selected from a face 
of tiie probe referred to as the "probe face plane." The user may create a transparent cut 
probe as described above to interface with an opaque ribbon section. Alternatively, the 
probe can be made opaque and the ribbon section can be made transparent if desired. 

The control points may be used to produce a plurality of line segments such as 
line segments 1722, 1724, and 1726, which collectively are referred to as polyline 1728 
wliich is like a polygon but may or may not be closed. Therefore, tlie line segments may 
form an open or closed line so that a single or multiple coolde planes may be produced. 
In a preferred embodiment, the area of data display, i.e. the coolde planes, is projected 
along a direction perpendicular to the probe fece plane and the data displayed may extend 
to an opposing face of the probe referred to as the "opposing probe face plane." After 
creating a ribbon section, the user may use the mouse controls or keyboard to select, 
move, drag, or grab the control points to edit the ribbon section 1710 in real time and 
display different data from the 3D volume data set along the cookie planes 1712 and 
1714. In addition to editing the ribbon section 1710 to view different data from the 3D 
volume data set within the probe, the entire ribbon section 1710 and probe may be 
simultaneously moved to a different position m order to view different data from the 3D 
volume data set outside the boundaries of the probe at its prior position. Active control 
point 1720 is preferably high-lighted or colored differently as compared to the other 
control points to indicate control point 1720 is in an active state for operations such as 
moving, deleting, or otherwise editing as discussed further below. Control points may 
be inserted or deleted after the initial ribbon section construction. The ribbon section 
geometry and orientatioii may also be saved for future work sessions. 

FIG. 18 depicts a block diagram of system 1810 for program modules in a 
presently preferred embodiment of the invention for producing ribbon sections at real 
time frame rates as discussed herein above. Probe module 1822 provides initial context 
for sketching polyline 1728. Probe module 1822 suppUes user activity data such as 
mouse clicks or keystrokes to cookie manager 1824. Thus, user activity data such as 
creation of control points, deletion of control points, moving of control points, moving 
of the entire probe, and the like are inserted into system 1810. 
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Cookie manager 1824 manages the user input data supplied by probe 1822. 
Cookie manager 1824 distributes the data, e.g. control points add, move, and delete as 
appropriate to both polyline module 1826 and cookie plane module 1828. For some 
types of data, such as inserting a control point, cookie manager 1824 receives data from 
polylme module 1826 and passes the data to cookie plane module 1828. 

Polyline module 1826 manages data related to polyline 1728 and the associated 
control points in conjunction with polymarker module 1832 and polystate module 1830. 
Polyline 1728 is mainly provided for visual reference. Polyline module 1826 in 
conjunction with polystate module 1830 manages the state of the control points. For 
instance, in the active state, a control point can be moved or deleted. The active control 

point, suchasactivecontrolpointmo is preferably highKghteA The controlpoints may 
be enlarged or decreased in size for easier viewing. Polymarker module 1832 provides 
visual context such as highlighting or varied colore, for the control pomts so that an 
operator knows which point is m the active state for moving, deleting, and otherwise 
editing. Polymarker module 1832 also may pro^dde text such as the location mdication 
adjacent active control point 1 720. 

Cookie plane module 1828 provides textured geometry, which may for mslance be 
lithography-related for geological data, to the surface of the cookie planes such as cookie 
planes 1712 and 1714. Cookie state module 1834 monitors the state of the cooMeplanes 
so fliat m the active state one or more cookie planes can be moved or otherwise edited 
whereas in flie inactive state no changes are made. 

During operationof system 1810, prt>be module 1822 may notifycookie manager 
1824 that an event has taken place, e.g., a marker deletion, i.e.. a control point deletion. 
Cookie manager 1824 then notifies polylme module 1826 which deletes the marker or 
control pomt and joins up the two surrounding lines into one line and notifies the 

polymarkermodule to removethedeleted marker firom the Ustofpolymarkersmaintamed 
by polymarker module 1832. Cookie manager 1824 and cookie plane module then 
convert (he two planes mto one. 

FIG. 19 iUustrates one embodiment of the invention for creating a three 
dimensional surfece representative of a physical phenomena described by a 3D volume 
data set such as, for example, a geological fault surfece. hi fact, this method has been 
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found to easily and quickly provide 3D volume visualization for rapid identification and 
interpretation of geological fault surfaces. However, other types of surfaces for other 
types of data could also be quickly described from use of the method of the present 
invention. Although the method is described in more detail below, in general the method 
5 is used to digitize control points, such as control points 1902, 1904, 1906, and 1908 along 
aprobeface 1910. Visual examination ofthe textured surface of probe face 1910pennits 
an operator to visually locate control points on a structure of interest such as a suspected 
fault line. Probe face 1910 is then moved, and a new set of control points are, digitized. 
Surface 1912 is then interpolated between the initial points and the new control points. 
1 0 The control points may be easily edited or moved to more accurately define the surface 
may be 1912 at a real time firame rate as discussed above. This process can be repeated 
until the surface interpretation is complete, at which time the surface 1912 can be saved. 

Thus, one embodiment of the invention disclosed by FIG. 1 9 provides a method 
of rapidly constructing a three-dimensional surface or fault found within a 3-D volume 
15 data set. In a preferred embodiment, 3-D probes, as discussed above, are used. The 
method involves constructing a plurality of spline curves, such as spline curve 1914 on 
the probe face 1910 which may be interactively constructed when visualizing the data 
displayed on probe face 1910. Spline curve 1914 is interpolated using an algorithm 
created in the same way as spline curves 1916, 1918, and 1920 are interpolated. Spline 
curves and v-curves 1922, 1924, and 1926, are produced to form a grid. The grid outlines 
the three dimensional sur&ce 1912. 

For construction ofthe initial spline curve, such as for instance spline curve 1914, 
the user digitizes control points, such as control points 1902, 1904, 1906 and 1908 on the 
probe face 1910. Markers are produced at these control points and spline curve 1914 is 
interpolated between the control points 1902, 1904, 1906 and 1908. Control points 1902, 
1904, 1906 and 1908 may be moved within probe face 1910, thereby interpolating a new 
spline curve 1914. Once the probe face 1910 is selected, the other probe surfaces are 
made transparent for ease of operation. Additionally, the selected probe face 1910 may 
be made opaque in order to view surface 1912 through the probe face 1910. 

The user then moves the probe face to 1910 to another position and selects new 
control points. The user may easily move back and forth between previously created 
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spline curves by selecting grid intersections such as intersection 1928 or 1930. As 
additional spline curves are created, the v-curves may also be smoothly and quickly 
inteipolated using another algorithm. The user preferably creates a plurality of spline 
curves in the same manner, and interpolation of surfece 1912 is immediately displayed 
at real time frame rates as discussed here above. 

The user may stop the probe and move one or more control points, such as control 
points 1902, 1904, 1906 and 1908 to adjust the position of the respective spline curve 
such as spline curve 1914. All other spline curve, such as spline curves 1916, 1918, and 
1920 remain the same while surface 1912 is smoothly interpolated between the current 
spline curve 1914 on the probe fece 1910 andthe prior spline curve 1920. The remainder 
of the surfece 1912 remains the same unless the user moves ±e prove face 1910 to 
another spline curve such as 1920 and proceeds to edit the same thereby reshapmg the 
surfece 1912 between the current spline curve 1920 on the probe face place 1910 and the 
prior splme curve 1918. 

Additional spline curves may be added between existing spline curves if desired. 
By selectmg grid intersections such as grid intersections 1928 or 1930, or by selecting 
control points on the probe face, such as control points 1902, 1904, 1096 and 1908, the 
user may quickly move the probe face 1910 back and forth as desked. Once a grid 
intersection is selected, the user may move therespective control points within the probe 
fece 1910 that will be displayed as indicated in HQ. 19, and suifece 1912 with its spline * 
curves and v-curves will follow interactively. The spline curves and the v-curves may 
or may not be displayed depending on preference of the user. Preferably, only one probe 
fece 1910 is displayed at a time for clarity. 

FIG. 20 discloses modular system 2000 which describes the software for 
performing the functions as described in connection with FIG. 19. In one embodiment, 
modular system 2000 incorporates many of the same modules as used in modular system 
1810 forproducingaribbon section. Thus, the design may permit a more general control 
that may effectively perform both ribbon section functions as weU as surface mapping 
functions. For instance, cookie module 2018, polymarker module 2012, and polyline 
module 2014 may be used in the manner discussed above relating to producing ribbon 
sections. Likewise, probe module 2010 therefore again performs the functions as 
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previously noted such as providing initial sketching splines such as spline 1914, which 
will also be observed to be similar to the polylines discussed hereinbefore. Probe module 
2010 supplies user activity data such as noiouse clicks and keystrokes to the various other 
modules. Cubic spline module 2016 relates to controlling functions for the editable 
suif ace or mapping feature such as surface 1912 shown in FIG. 1 9. Cubic spline module 
2016, spline manager module 2020, spline curve module 2022, and spline surface module 
2024 ai-e functions that are different from those used only in creating ribbon sections. 
Spline manager 2020 performs a number of different functions such as creating or 
deleting surfaces, changing mode functions from creating a surface to editing a surface, 
as well as read and write functions. For instance, spline manager 2020 may read and 
write to a surface attribute file regarding attributes such as colors, showing the grid, 
marker colors, and the lilce. Spline curve module 2022 keeps track of die splhie curves, 
keeps track of which splme curve may be m the editing or creatuig state and keeps track 
of changes made to the spline curves. Splme surface module 2024 acts on the changes 
made for interpolating the resulting surface changes. 

By using the system and method of the present invention, geologists and 
geophysicists can more quickly and accurately visualize and interpret 3D seismic data. 
This sharply reduces 3D seismic project cycle time, boosts production from existing 
fields, and finds more reserves. 

Conclusion 

While various embodiments of the present invention have been described above, 
it should be understood that ttiey have been presented by way of example only, and not 
limitation. Thus, the breadth and scope of the present invention shovdd not be limited by 
any of the above-described exemplary embodiments, and should be defined in accordance 
v^th the spirit of the mvention including the following claims and their equivalents. 
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Claims 

What Is Claimed Is: 

1. A method for imaging a three-dimensional data volume, said three- 
dimensional data volume comprising apluralitj^ of voxels, each voxel comprising a three- 
dimensional location and a dataword, said dataword being representative of a physiced 
phenomena, said method comprising; 

creating at least one tluee-dimensional sampling probe, wherein said at 
least one three-dimensional sampling probe is the same size or a subset of said 
three-dimensional data volume, said at least one three-dimensional sampling 
probe having a probe face in a probe face plane and an opposing probe face in an 
opposing probe face plane; 

producing a plurality of control points in said probe face plane, said 
plurality of control points defining one or more lines on said probe face plane; 

extendmg a ribbon section firom said probe face plane toward said 
opposing probe face plane, one edge of said ribbon section being formed by said 
one or more lines; and 

selectively imaging datawords representative of said physical phenomena 
at three-dimensional locations which intersect said ribbon section and said three- 
dimensional sampUng probe. 

2. The method of Claim 1, further comprising: 

editing said plurality of control points in said face plane to thereby 
redefine said one or more lines, and 

extending another redefined ribbon section fi-om said probe face plane 
toward said opposing probe face plane. 

3 . The metiiod of Claim 2, wherein said step of editing further comprises: 
deleting one or more of said plurality of control points. 

4. The method of Claim 2, wherein said step of editing further comprises: 
changing a location of one or more of said plurality of control points. 
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5. The method of Claim 2, wherein said step of editing further comprises: 
adding one or more control points to said plurality of control points. 

6. The method of Claim 1, wherein said ribbon section is perpendiculai* to 
said probe face plane. 

7. The method of Claim 1, wherein said ribbon section extends from said 
probe fiice plane to said opposing probe face plane. 

8. The method of Claim 1, wherein said one or more lines comprise a 
plurality of straight lines. 

9. The method of Claim l,whereinsaid one or more lines form a closed line. 

10. The method of Claim 1, wherein said ribbon section is comprised of a 
plurality of planes. 

11. The method of Claim 1, wherein said three-dimensional probe has a 
plurality of side faces perpendicular to said probe face plane, said ribbon section being 
unparallel with respect to each of said plurality of side feces. 

12. A program storage device readable by a machine, embodying a program 
of instructions executable by machine to perform method steps for imaging a three- 
dimensional data volume, said three-<limensional data volume comprising a plurality of 
voxels, each voxel comprising a three-dimensional location and a dataword, said 
dataword being representative of a physical phenomena, said method comprising; 

displaying a plane within said three-dimensional data volume; 
producing a plurality of control points in said plane, said plurality of 
control points defining one or more lines on said plane; 
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extending a ribbon section from said plane, one edge of said ribbon 
section being formed by said one or more lines; and 

selectively imaging said datawords representative of said physical 
phenomena at three-dimensional locations which intersect said ribbon sectionand 
said three-dimensional data volume. 

1 3 . Hie method of Claim 12, further comprising: 

editing said pluraUty of control points in said plane to thereby redefine 
said one or more lines, and 

extending another redefined ribbon section firom said plane toward an 
opposing plane. 

14. The method of Claim 13, wherein said step of editing further comprises: 
deleting one or more of said plurality of control points. 

15. The method of Claim 13, wherein said step of editing fixrther comprises: 
changing a location of one or more of said plurality of control points. 

16. ThemethodofClaiml3,whereinsaidstepofeditingfiarthercomprises: 
adding one or more control points to said pluraUty of control points. 

17. ThemethodofClaiml2,wbereinsaidribbonsectionispeipendicularto 
said plane. 

18. ThemethodofClaim 12, wherein said ribbon section extends from said 
probe plane to an opposing plane. 

19. The method of Claim 12, wherein said one or more lines comprise a 
plurality of straight lines. 
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The method of Claim 12, wherein said one or more lines fonn a closed 



2 1 . The method of Claim 1 2, wherein said ribbon section is comprised of a 
plm^ality of planes. 

22. The method of Claim 12, wherein plane is a face plane of a three- 
dimensional probe which also has a plurality of side faces perpendicular to said face 
plane, said ribbon section being oriented so as not to be parallel to one or more of said 
plurality of side faces. 



23. A method for imaging a three-dimensional data volume, said three- 
dimensional data volume comprising a plurality of voxels, each voxel comprising athree- 
dimensional location and a dataword, said method comprising: 

positioning a face of a probe at a &st position within said three- 
dimensional data volume; 

forming a &st set of control points on said face of said probe for tracking 
aphysicalphenomenadescribedby said threeslimensional data volume, said first 
set of control points defining a first spline curve- 
moving said face of said probe to a second position within said three- 
dimensional volxmie; 

forming a second set of control points on said fece of said probe for 
tracking said physical phenomena, said second set of control pomts defining a 
second spline curve; and 

interpolating between said first spline curve and said second spline curve 
to define a three dimensional surface representative of said physical phenomena. 
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24. The method of Claim 23; further comprising: 

displaying the surface representative of said physical phenomena, said 
surface intersecting said first set of control points and said second set of control 
points. 

25. The method ofCIaim 23, fturther comprising: 

interpolating between said first set of control points to define said first 
spline curve and interpolating between said second set of contiol points to defme said 
second spline curve, at least one of said first spline curve and said second spline curve 
being cimdlinear 

26. The metiiod of Claim 23, further comprising: 

moving said fece of said probe to a third position within said three- 
dimensional volume; 

forming a third set of control points on said face of said probe for tracking 
said physical phenomena, said third set of control points defining a third spline 
curve; and 

interpolating between said first spline curve, said second spline, and said 
third spline curve for enlarging said surface. 

27. The method of Claim 23, fijrther comprising: 

editing at least one of said first set of control points and said second set 
of control points. 



28 . The melhod of Claim 23, fiirther comprising: 

displayingv-curves which interconnect betweenrespectivecontrolpoints 
at said first position of said probe and said second position of said probe. 
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29. The method of Claim 28, further comprising: 

displaying said spline curves and said v-curves, said spline' curves and 
said v-curves forming a grid representative of said physical phenomena, said grid 
having a plurality of intersections between said spline curves and said v-curves. 

30. The method of Claim 29, further comprising: 

selecting one of said plurality of intersections, and moving said 
intersection to thereby edit said grid. 



3 1 . The mefliod of Claim 29, further comprising: 

selecting one of said plurality of intersections to thereby reposition said 
face to pass through said mtersection. 



32. The method of Claim 3 1 , further comprismg: 

selecting one of said second set of control points to thereby reposition said 
face to pass through said second set of control points. 

33 . A program storage device readable by a machine, embodying a program 
of instructions executable by machine to perform method steps for imaging a three- 
dimensional data volume, said three-dimensional data volume comprising a plurality of 
voxels, each voxel comprising a three-dimensional location and a dataword, said method 
comprising: 

positioning a plane at a plurality of plane positions within said three- 
dimensional data volume; 

forming a set of control points at each of said plurality of plane positions 
such that each of said set of control points defines a related spline curve; and 

interpolating between each of said spline curves to form a surface 
representative of a said physical phenomena described by said three-dimensional 
data volume. 

34. The method of Claim 33, further comprising: 
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displaying said surface representative of said physical phenomena, said 
surface intersecting each of said set of control points. 

35. The method of Claim 33, further comprising: 

interpolating between each of said set of control points to define 
said related spline curves, at least one of said related spline curves being 
curvilinear. 



36. The method of Claim 33, further comprising: 
editing one or more of said control points. 

37. The method of Claim 33, further comprising: 

displaying v-curves which interconnect between respective control points 
at said plui-ality of plane positions. 

38. The method of Claim 37, further comprismg: 

displaying said spline curves and said v-curves to form a grid 
representative of said physical phenomena, said grid having a plurality of 
intersections between said spline curves and said v-curves. 

39. The method of Claim 38, further comprising: 

selecting one of said pluraUty of intereections, and moving said 
intersection. 



40. The method of Claim 39, fiirther comprising: 

selecting one of said plurality of intersections to thereby reposition said 
plane to pass througli said intersectioa 

41. The method of Claim 40, further comprising: 

selecting one of said control points to thereby reposition said plane to pass 
through said selected one of said control pomts. 



-48- 



wo 02/37422 



PCT/USOO/29835 



AMENDED CLAIMS 
[received by the International Bureau on 25 March 2002 (25.03.02); 
original claims 1-41 replaced by new claims M3 (8 pages)] 

1. A method for imaging a three-dimensional data volume, said three- 
dimensional data volume comprising a plurality of voxels, each voxel comprising a three- 
5 dimensional location and a dataword, said dataword being representative of a physical 
phenomena, said method comprising; 

creating at least one three-dimensional sampling probe, wherein said 
three-dimensional sampling probe is the same size or a subset of said three- 
dimensional data volume, said three-dimensional sampling probe having a probe 
face in a probe face plane and an opposing probe face in an opposmg probe face 
plane; 

producing a plurality of control points in said probe face plane, said 
plurality of control points defining one or more lines on said probe face plane; 
extending a ribbon section from said probe face plane toward said 
^ ^ opposing probe face plane, one edge of said ribbon section being formed by said 

one or more lines; and 

selectively imaging datawords representative of said physical phenomena 
at three-dimensional locations whicti intersect said ribbon section and said three- 
dimensional sampling probe. 

20 

2. The method of Claim 1 , further comprising: 

editing said plurality of control points in said probe face plane to thereby 
redefine said one or more lines, and 

extending another redefined ribbon section from said probe face plane 
25 toward said opposing probe face plane. 

3. The method of Claim 2, wherein said step of editing further comprises: 
deleting one or more of said plurality of control points. 



30 



4. The method of Claim 2, wherein said step of editing further comprises: 
changing a location of one or more of said plurality of control points. 
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5. The method of Claim 2, wherein said step of editing further comprises: 
adding one or more control points to said plurality of control points. 

6. The method of Claim 1, wherein said ribbon section is perpendicular to 
said probe face plane. 

7. The method of Claim 1 , wherein said ribbon section extends from said 
probe face plane to said opposing probe face plane. 

8. The method of Claim 1, wherein said one or more lines comprise a 
plurality of straight lines. 

9. The method of Claim 1, wherein said one or more lines form a closed line. 

10. The method of Claim 1, wherein said ribbon section is comprised of a 
plurality of planes. 

11. The method of Claim 1, wherein said three-dimensional probe has a 
pluraUty of side faces perpendicular to said probe face plane, said ribbon section being 
unparallel with respect to each of said plurality of side faces. 

12. A program storage device readable by a machine, embodying a program 
of instructions executable by machine to perform method steps for imaging a three- 
dimensional data volume, said three-dimensional data volume comprising a pluraUty of 
voxels, each voxel comprising a three-dimensional location and a dataword. said 
dataword being representative of a physical phenomena, said method comprising; 

displaying a plane within said three-dimensional data volume; 
producing a plurality of control points in said plane, said pluraUty of 
control points defining one or more lines on said plane; 
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extending a ribbon section from said plane, one edge of said ribbon 
section being formed by said one or more lines; and 

selectively imaging said datawords representative of said physical 
phenomena at three-dimensional locations which intersect said ribbon section and 
5 said three-dimensional data volume. 



13. The method of Claim 12, further comprising: 

editing said plurality of control points in said plane to thereby redefine 
said one or more lines, and 

extending another redefined ribbon section from said plane toward an 
opposing plane. 



14. The method of Claim 13, wherein said step of editing further comprises: 
deleting one or more of said plurality of control points. 

15 

15. The method of Claim 13, wherein said step of editing further comprises: 
changing a location of one or more of said plurality of control points. 

16. The method of Claim 13, wherein said step of editing further con:q)rises: 
20 adding one or more control points to said plurality of control points. 

17. The method of Claim 1 2, wherein said ribbon section is perpendicular to 
said plane. 

^5 18. The method of Claim 12, wherein said ribbon section extends fi-om said 

plane to an opposing plane. 



19. The method of Claim 12, wherein said one or more lines comprise a 
plurality of straight lines. 

30 
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20. The method of Claim 12, wherein said one or more lines form a closed 

line. 



21. The method of Claim 12, wherein said ribbon section is comprised of a 
plurality of planes. 

22. The method of Claim 12, wherein said plane is a probe face plane of a 
three-dimensional probe which has a plurality of side faces peipendicular to said probe 
face plane, said ribbon section being oriented so as not to be parallel to at least one of 
said plurality of side faces. 

23. A method for imaging a three-dimensional data volume, said three- 
dimensional data volume comprising a plurality of voxels, each voxel comprising a three- 
dimensional location and a dataword, said method comprising: 

positioning a face of a probe at a first position within said three- 
dimensional data volume; 

forming a first set of control points on said face of said probe for tracking 
aphysical phenomena describedbysaid three-dimensional data volume, saidfirst 
set of control points defining a first spline curve; 

moving said face of said probe to a second position within said tiuee- 
dimensional volume; 

forming a second set of conteol points on said face of said probe for 
tracking said physical phenomena, said second set of control points defining a 
second spline carve; and 

interpolating between said first spline curve and said second spline curve 
to define a three-dimensional surface representative of said physical phenomena. 
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24. The method of Claim 23, further comprising: 

displaying the surface representative of said physical phenomena, said 
surface intersecting said first set of control points and said second set of control 
points. 



25. The method of Claim 23, further comprising: 

interpolating between said first set of control points to define said first 
spline curve and interpolating between said second set of control points to define said 
second spline curve, at least one of said first spline curve and said second spline curve 



being curvilinear. 



26. The method of Claim 23 , further comprising: 

movmg said face of said probe to a tiiird position within said three- 
dimensional volume; 

forming a third set of control points on said face of said probe for tracking 
said physical phenomena, said third set of control points defining a third spline 
curve; and 

interpolating between saidfiirst spline curve, said second spline curve, and 
said third spline curve for enlarging said surface. 

27. The method of Claim 23, further comprising: 

editing at least one of said first set of control points and said second set 
of control points. 

28. The method of Claim 23, further comprising: 

forming a plurality of v-curves which interconnect between respective 
control points at said first position of said probe and said second position of said 
probe. 
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29. The method of Claim 28, further comprising: 

displaying said spline curves and said v-curves, said spline curves and 
said v-curves forming a grid representative of said physical phenomena, said grid 
having a plurality of intersections between said spline curves and said v-curves. 

5 

30. The method of Claim 29, further comprising: 

selecting one of said plurality of intersections, and moving said 
intersection to thereby edit said grid. 

10 3 1 . The method of Claim 29, further comprising: 

selecting one of said plurality of intersections to thereby reposition said 
face to pass through said intersection. 

32. The method of Claim 29, further comprising: 

selecting one of said first set of control points and said second set of 
control points to thereby reposition said face to pass through one of said first set 
of control points and said second set of control points. 



20 
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33. The method of Claim 23, further comprising: 

forming a third set of control points on said face of said probe at said first 
position, said third set of control points defining a third spline curve; 

forming a fourth set of control points on said face of said probe at said 
second position, said fourth set of control points defining a fourth spline curve; 
and 

interpolating between said third spline curve and said fourth spline curve 
to define another three-dimensional surface representative of another physical 
phenomena described by said three^mensional data volume, said three- 
dimensional surface and said another three-dimensional surface being defined 
substantially at (he same time. 



30 
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34. A program storage device readable by a machine, embodying a program 
of instructions executable by machine to perform method steps for imaging a three- 
dimensional data volume, said three-dimensional data volume comprising a plurality of 
voxels, each voxel comprising a three-dimensional location and a dataword. said method 
comprising: 

positioning a plane at a plurality of plane positions u^ithin said three- 
dimensional data volume; 

forming a set of control points at each of said plurality of plane positions 
such that each of said set of control points defines a related spline curve; and 

interpolating between each of said spline curves to fomi a surface 
representative of a physical phenomena described by said three-dimensional data 
volume. 



35. The method of Claim 34, further comprising: 

displaying said surface representative of said physical phenomena, said 
surface intersecting each of said set of control points. 

36. The method of Claim 34, further comprising: 

interpolating betv^^een each of said set of control points to define said 
related spline curves, at least one of said related spline curves being curvilinear. 

37. The method of Claim 34, further comprising: 
editing one or more of said control points. 

38. The method of Claim 34, further comprising: 

forming a plurality of v-curves which interconnect between respective 
control points at said plurality of plane positions. 
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39. The method of Claim 38, further comprising: 

displaying said spline curves and said v-curves to fonn a grid 
representative of said physical phenomena, said grid having a plurality of 
intersections between said spline curves and said v-curves. 

5 

40. The method of Claim 39, further comprising: 

selecting one of said pluraUty of intersections, and moving said 
intersection to thereby edit said grid. 

^® "^1- The method of Claim 39, further comprising: 

selecting one of said plurality of intersections to thereby reposition said 
plane to pass tiirough said intersection. 



15 
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42. The metiiod of Claim 39. fiirther comprising: 

selecting one of said sets of control points to thereby reposition said pi; 
to pass through said one of said sets of control points. 



ane 



43. The method of Claim 34. further comprising: 

forming another set of control points at each of said plurality of plane 
positions, such that each of said anotiier set of control points defines another 
related spline curve; and 

interpolating between each of said anotiier spline curves to form anotiier 
surface representative of anotiier physical phenomena described by said fliree- 
dimensional data volume, said surface and said anotiier surface being formed 
25 substantially at tiie same time. 
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